Interaction of discrete nonlinear Schrodinger solitons with a linear lattice impurity 
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The interaction of moving discrete solitons with a linear Gaussian defect is investigated. Solitons 
with profiles varying from hyperbolic secant to exponentially localized are considered such that the 
mobility of soliton is maintained; the condition for which is obtained. Studies on scattering of the 
soliton by an attractive defect potential reveal the existence of total reflection and transmission 
windows which become very narrow with increasing initial soliton amplitude. Transmission regions 
disappear beyond the small-amplitude limit. The regions of complete reflection and partial capture 
correspond to the windows of the existence and nonexistence of solution of the stationary problem. 
Interaction of the discrete soliton with a barrier potential is also investigated. The critical amplitude 
of the defect at which splitting of the soliton into two parts occurs was estimated from a balance 
equation. The results were confirmed through direct numerical integration of the dynamical equation 
showing very good agreement with the analytical prediction. 



I. INTRODUCTION 

Since the first theoretical prediction [l| and experimen- 
tal observation of discrete sohtons in photonic crystals 
there has been a lot of activity related to the gen- 
eration, stability and dynamical properties of locahzcd 
modes in discrete nonlinear systems (see reviews 0] and 
references therein) . Another physical system where a dis- 
crete setting naturally appears is Bose-Einstein conden- 
sates (BECs) embedded in ID, 2D or 3D optical lattices 
[3] • For some specific parameters of the system (relatively 
high amplitude of the optical lattice) the Gross-Pitaevskii 
equation which governs the dynamics of the continuum 
wave function can be approximated by a discrete nonlin- 
ear Schrodinger (DNLS) equation @. The existence of 
discrete solitons has been demonstrated in a wide range 
of other physical systems such as, for example, atomic 
chains with on-site cubic nonlincarities biophysical 
systems Q, and in arrays of coupled nonlinear optical 
wave guides 0]. Very recently discrete gap solitons were 
also observed experimentally in a saturable media with 
alternate spacing of waveguide arrays Q . 

In this work, we will focus on the propagation of dis- 
crete ID solitons through a nonlinear lattice in the pres- 
ence of a localized linear defect. The classical process of 
scattering of an incident wave (particle) on a localized 
defect in the form of a barrier/well is well known and in- 
vestigated. It basically consists of three different regimes; 
reflection, transmission and capture depending on the ra- 
tio between kinetic energy of the incoming wave/particle 
and the height of the defect. However treating the same 
process as a quantum mechanical problem results in in- 
teresting phenomena. Interaction of a discrete soliton 
with a localized defect in the intcgrable Ablowitz-Ladik 
model and DNLS equation was considered theoretically 
in j9l-[l3|. The first experimental attempt to describe 
the discrete soliton interaction with a localized defect in 
arrays of coupled waveguides revealed significant differ- 
ences in the linear and nonlinear regimes of the scatter- 
ing process [l3|- A great deal of attention was paid to 



the study of soliton interaction with a (5-likc linear defect 
in the continuum nonlinear Schrodinger (NLS) equation 
[isl [l6j . Recently this problem regained strong inter- 
est due to the possibility of using the interaction of the 
matter-wave solitons with a defect in BEG to investi- 
gate different quantum effects as well as to study the 
interaction of BEG with surfaces [13 ■ Matter- wave soli- 
ton interaction with extended defect was investigated in 
Ref . [l8l - l22j where the authors discussed the phenomenon 
of resonant transmission and trapping of bright solitons 
in the presence of a quantum well. Interaction of gap soli- 
tons with a localized defect in an optical lattice [23, [IS] 
and trapping of discrete solitons by linear and nonlin- 
ear defects in different discrete models (2B - [29| have also 
been considered. Scattering of a single matter-wave soli- 
ton and two-soliton molecule in a potential well was also 
recently considered in (soj . 

The article is organized as follows: the model is intro- 
duced and discussed in section pH where we have also 
obtained the frequency condition for the mobility of the 
soliton, allowing us to consider both hyperbolic secant 
and exponentially localized solitons which are not de- 
stroyed or pinned to the lattice. Section (|III[) discuses 
the scattering process for an attractive defect potential, 
which is followed by a defect mode analysis in section 
PV[) . Finally, splitting of discrete soliton by a repulsive 
defect is considered in section dVl). 



II. MODEL 

One of the basic lattice models which appeared in the 
last decades in various contexts of physics and biology 
and has been intensively studied (review [^) is the DNLS 
equation, which has the form 

^~tr ^ ^2^" + 2cr I V'ri 1^ (1) 

where ipn = V'n(0 is the time dependent amplitude of the 
individual waveguide modes when an array of waveguides 
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is considered, and is the amplitude of the individual os- 
cillator at site n when a lattice of coupled anharmonic 
oscillators is considered, A2'0n = ipn+i +'0n-i — de- 
scribes the coupling between neighboring waveguides (or 
lattice sites), cr = ±1 is the nonlinear coefficient which 
accounts for the optical Kerr effect (or interparticlc in- 
teraction in BEC), and the linear potential Vn takes into 
account the presence of a localized impurity in the lattice. 

In the absence of the impurity, Vn = 0, the simplest 
class of discrete breathers can be obtained using the 
steady-state ansatz ipnit) = e'"*u„, where lj corresponds 
to the frequency of the solution. Then the stationary am- 
plitudes Un, n = 0, ±1, ±2... satisfy the lattice equation: 

A2U,i + LdUn +2a \un P u„ = (2) 

with the condition for spatial localization u„ — > as rt — s- 
±00. 

The dispersion relation, which describes the propaga- 
tion of linear waves, and can be obtained considering 
plane wave solution of the form u„(t) oc e**"'", is well- 
known: 

w = 4sin2(fc/2), (3) 

and the frequency is bounded in the interval < a; < 4 
(linear band). Above and below the linear band there are 
two semi-infinite gaps where in the presence of nonlinear- 
ity one can excite nonlinear localized modes. Such modes 
exist due to the balance between dispersion and nonlin- 
ear effects. It can be shown that the solutions in these 
semi-infinite gaps are connected through the following re- 
lation: if the discrete function is the solution from 
the semi-infinite gap below the linear band with a; < 
and (7=1 then the function {(— 1)"m„} is the solution of 
Eq. for (T = — 1 and lu > 4. Therefore in the following, 
we will work only with localized solutions in one of the 
semi-infinite gap, namely, below the linear band a; < 0, 
where a = 1. 

Inside the semi-infinite gap one can consider two limit- 
ing cases allowing for an analytic description of the fun- 
damental discrete soliton. In the vicinity of the linear 
band edge the problem can be mapped to the contin- 
uum NLS equation with onsite nonlincarity or can be ob- 
tained as the tight-binding approximation of the Gross- 
Pitaevskii equation in the presence of deep optical lattice. 
Essentially, that leads to the sech-type approximation of 
the small-amplitude (envelope) discrete soliton: 

Un « A sech(^n), (4) 

where A and 1/A are the amplitude and the character- 
istic width of the discrete soliton, respectively. For this 
solution, the total power P and A has a linear depen- 
dence: 

n 

In the opposite limit, when the soliton frequency lies 
far from the linear band with ^ ^ 1, the solution be- 
comes highly localized (intrinsic localized mode), and an 
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FIG. 1. (Color online) (a) Power dependence of the defect- free 
fundamental soliton on the central peak amplitude. Discrete 
soliton profile for (b) uj = —0.05 and (c) uj — —1. The red 
boxes, black solid line and blue dashed line correspond to nu- 
merical, sech-profile Q and exponential-profile ((6| solutions, 
respectively. 



exponential localization can be used to approximate the 
profile of the solution: 

u„«Ae-'^l"l. (6) 

In Figllja) , numerically found dependence of the total 
power P of the fundamental discrete soliton on its am- 
plitude is shown. In panels (b) and (c), the profiles of 
the fundamental discrete solitons calculated numerically 
with the two approximations ^ and (jS]) are compared. 
As we can see, in the vicinity of the linear band the small- 
amplitude approximation can be used to fit the numerical 
solution (see Fig[T](b)) while for large |a;|, where A>1, 
the numerical solution is well approximated by the expo- 
nential function (shown in Figll](c)). 

As we are interested in the dynamical interaction of 
the discrete soliton with a linear defect, it is necessary 
to have a stable mobile solution to observe its long-time 
dynamics. Here we discuss the limits at which the prob- 
lem of interaction of discrete soliton with defect can be 
considered which is related to the mobility of the soli- 
tons. It is known that the mobility of discrete solitons is 
restricted by the presence of the "Peierls-Nabarro" (FN) 
barrier which appears due to the non-integrability of the 
DNLS equation jlllll. Effectively, the FN barrier can 
be seen as the energy difference between on-site and inter- 
site fundamental modes with equal powers. Larger this 
difference stronger is the "friction" between soliton and 
the lattice, which leads to dissipation of the soliton by 
radiation during propagation. At a critical amplitude of 
the soliton the FN barrier becomes strong enough to stop 
the soliton and pin it to the lattice. The mobility of dis- 
crete solitons in terms of BEC in optical lattice under 
realistic experimental conditions was considered in fssj . 
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with frequency Ws , placed far from the defect at the posi- 
tion Hs, and moving towards the impurity with a velocity 
V. In order to estimate the shift in the soliton frequency 
due to the applied velocity, we will use a Taylor expan- 
sion of the dispersion relation ([3]) in the vicinity of the 
corresponding edge of the linear band modes fcp 

cj(fc) = uj{ko) + VB\k=ko{k~ko) + ^D\k^ko{k~ko)^, (8) 

where the Bloch velocity vb and diffraction coefficient 
D (equivalent to the inverse effective mass in solid-state 
physics) are defined as follows: 



2sin(A:), 



2cos(fc). (9) 
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Notice that Eq. ([8]) is valid only near the stationary 
points fco (the bottom edge of the linear band mode corre- 
sponds to fco = ±27rj and the top one to fco = i(2j + l)7r 
where j = 0,1,2,...), where the Bloch velocity can be 
approximated by vsik) ~ D{k — fco). Therefore, by ap- 
plying the initial velocity v = vb the frequency of the 
moving soliton ujs{v) can be estimated according to the 
following relations: 



FIG. 2. (Color online) (a) Time dependence of the maxi- 
mum amplitude maa;(|«„p) and (b) Trajectory for different 
frequencies uj = —0.05; —0.3; —0.4; —0.45; —0.5 of the soliton 
with an initial velocity v — 0.1. 



In order to find the frequency limit on mobility, we 
considered solitons with amplitude ranging from low to 
high, and studied their time dynamics for a fixed initial 
velocity. As shown in Figl2l three different regimes can 
be observed. Small amplitude solitons propagate freely 
maintaining their shape and without feeling the PN bar- 
rier (for Lu = —0.05). Whereas, high-amplitude solitons 
are pinned to the lattice after some distance of propaga- 
tion (w = —0.45 and —0.5). An intermediate regime is 
also observed, where the solitons are still able to propa- 
gate slowing down their velocities (w = —0.3 and —0.4)). 
This implies that, considering soliton solutions in the fre- 
quency range —0.4 < w < 0, allows one to study the 
interaction properties of both small-amplitude and high- 
amplitude solitons without the soliton being dissipated 
or pinned. 



III. SCATTERING PROCESS ON A 
POTENTIAL WELL, Vo > 

We begin by considering the interaction of the discrete 
soliton with an attractive defect (potential well) . As ini- 
tial soliton we consider the fundamental discrete soliton 
of the form 



r 



(10) 



where we have used the fact that in the vicinity of the 
edge of the linear band, fc = fco, it follows from Eq. © 
that VB{k = fco) = and D{k = fco) = 2. 

The localized defect is taken in the form of a Gaussian 
function 



y„ = yoexp[-(n-no)V(2d')] 



(11) 



where Vo < (Vb > 0) and d arc the amplitude and the 
width of the barrier (well), respectively, and tiq stands 
for the position of the impurity. In the following we will 
always consider the defect placed at the origin, hq = 0. 
A Gaussian profile of the defect can be easily introduced 
by external beams and can also reduce strong radiation 
processes [12, [13 which appear during interaction of the 
soliton with boundaries of the defect with sharp profiles. 

To analyze the soliton-defect interaction, we follow 
the dynamics of the reflection, capture and transmis- 
sion [RCT) coefficients which are defined in the following 
way: 



R 



C = 



T : 



n— — oo 
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p E \^r.\' 



P 



(12) 
(13) 
(14) 



-ivn/2-\-iuJs t 



(7) 



with the condition R + C + T = 1. Here the normaliza- 
tion constant P corresponds to the power of the initial 
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FIG. 4. Density plot for reflection (R), transmission (T) and 
capture (C) coefRcients as a function of (Vo, d) calculated as a 
solution of the dynamical equation with initial soliton taken 
at ljJs = —0.05 with initial velocity v = 0.1. 
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FIG. 5. Typical evolution dynamics for reflection (a), trans- 
mission (b) and capture (c). Initial and flnal soliton profiles 
are shown in top and bottom rows respectively. 



FIG. 3. (Color online) Dependence of the reflection (dashed 
black line), transmission (solid red line) and capture (blue 
dotted line) coefficients of the fundamental discrete soliton 
on the strength of the impurity, with the width d = 2, taken 
at (a) LJs = -0.05 (b) ujs = -0.1 (c) lu^ = -0.2 and (d) 
ujs = —0.3. The initial soliton velocity is u = 0.1 and is 
centered at Us = —50. 



stationary soliton. In order to determine the capture co- 
efRcient we take A equal to twice the width of the initial 
discrete soliton. 

The results for the RCT coefficients are obtained by 
the numerical integration of Eq. ([1]). A far placed soli- 
ton with a constant velocity is impinged on to the de- 
fect and is allowed to interact with the defect for a suf- 
ficiently long time. The RCT coefRcients arc calculated 
using Eq. p2p - (fTil) and this process is repeated for differ- 
ent potential depths Vq G [0,3]. The results are plotted 
in FiglSl The panels correspond to different frequencies 
of the initial soliton, such that the shape of the soliton 
changes from hyperbolic secant (low amplitude) to ex- 
ponential (high amplitude) but still in the regime where 
the soliton has mobility and is not affected strongly by 
the PN barrier. The results show that for some values of 
the depth of the defect Vq there are regions of complete 



reflection of the soliton by the potential well and also re- 
gions where transmission, reflection and partial capture 
occur. One of the interesting features of this result is 
the presence of sharp transmission lines where T « 1 
(red line in Fig|3]). The dependence of RCT coefficients 
as functions of Vq and d for ujs = —0.05 is plotted in 
FigHl As one can see, with increasing d and Vq, the re- 
gions where the soliton is reflected become narrower and 
windows of transmission become wider. The dynami- 
cal evolution of the small-amplitude soliton for the three 
different regimes; complete reflection, transmission and 
capture are presented in the middle row of FiglS] where 
the top and bottom rows correspond to the initial and 
final profiles of the soliton respectively. 

Similar diagrams for RCT coefficients on the inter- 
action of the gap soliton with Gaussian defect in the 
continuous Gross-Pitaevskii equation was published re- 
cently [2^. While generally the RCT diagrams from 
Fig|4] shows similar structure with Fig. 2 in [23| the be- 
havior of the transmission, reflection and capture curves 
for stronger defect are slightly different which can be ex- 
plained by different definition of the defects and also pres- 
ence of the optical lattice. 

The different frequencies in Fig|3] correspond to in- 
creasing amplitudes; A — 0.224 in (a), A = 0.318 in (b), 
A = 0.452 in (c) to ^ = 0.558 in (d). In all these cases the 
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IV. STATIONARY DEFECT MODE ANALYSIS 
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FIG. 6. (Color online) (a) Time dynamics of RCT coefficients 
and (b) corresponding soliton dynamics for uja — —0.2 and 
Vo = 1.2. 



FIG. 7. (Color online) Existence curve for the defect modes 
with different symmetries as a function of the strength of the 
potential for (a) ljs = —0.05 and (b) ujs = —0.3 for d = 2. 



soliton is able to propagate without being destroyed or 
pinned to the lattice (see Figl2|). One of the noticeable 
feature of the interaction of the high-amplitude soliton 
with the potential well is the reduction of the width of 
the total reflection regions. By increasing the amplitude 
of the soliton one also observes that the transmission re- 
gions disappear [see Figl3Jc), (d)]. This can be explained 
by the fact that in the high-amplitude limit the defect 
modes become highly unstable which prevents the dis- 
crete soliton from transforming into a stable defect mode 
and pass through the potential well. In order to show this 
we considered initial soliton at uig = —0.2 and Vb = 1.2, 
where the total reflection has step-like decrease and total 
transmission of the discrete soliton through the defect is 
expected. The dynamical picture of the RCT coefficients 
as well as evolution of the soliton profile are shown in 
Figini First, when the soliton enters the defect region it 
becomes totally captured as can be seen from the dynam- 
ics of the caption coefficient C{i) at t = 200 in FiglBI^a). 
It is visible that during the passage through the defect 
the soliton transforms into the on-site defect mode with 
three peaks (also in agreement with the stationary de- 
fect modes analysis). However, the soliton is not able to 
transform into a stable mode of the defect and instead of 
transmission, the soliton becomes trapped by the defect, 
continuously radiating the extra energy and transforming 
itself into a one peak stable solution. This kind of un- 
stable dynamics also explains the observed behavior of 
RCT coefficients between reflection regions in Fig|31Jc), 
(d) where strong radiation into the transmission and re- 
flection regions increases corresponding coefficients. 
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FIG. 8. (Color online) Profiles of impurity modes at LOa = 
—0.05 (black solid line, circles) and u)s ~ —0.3 (red dashed 
line, squares) with d = 2. In panels (B)-(D) the profiles of 
the defect modes at Vo = 0.39; 1.16; 2.1 for lj — —0.05 and 
Vb — 0.64; 1.5; 2.53 for cj = —0.3 with the same power as 
fundamental soliton P = Pa presented in (A) are shown. 

The existence of the regions of total reflection, total 
transmission and partial capture of the discrete soliton 
by the defect can also be understood by analyzing the 
stationary modes of the corresponding nonlinear system 
[2^. The stationary problem ^ in the presence of the 
linear defect now becomes 

A2U„ + LdgUn + VnUn + 2 \ Un ^ = 0. (15) 

Fixing the frequency of the soliton, lo^ , the width of the 
potential well, d, and changing the amplitude of the po- 
tential Vb, we found the existence of defect modes of 
two different types which are dynamically stable, namely, 
OS (on-site symmetric) and OA (on-site antisymmetric). 
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FIG. 9. (Color online) Dependence of the capture (blue dotted 
line) coefficient on the amplitude of the impurity in compar- 
ison with existence curves for the SO and AO defect modes 
normalized to the power of the fundamental soliton (P/Pa). 



The dependence of power of these defect modes with the 
strength of the potential with a fixed width (d = 2) for 
frequencies lo^ ~ —0.05 and —0.3 is presented in FiglT] 
The power of the dcfcct-frcc soliton, Pa, is shown as the 
horizontal line crossing the branches of the existence of 
defect modes at the matching points where Pa = Pj, 
j =B,C,D. In FiglSl the profiles of the fundamental soli- 
ton (A) and the first three defect modes at the match- 
ing points (B)-(D) are shown for both small and high- 
amplitude soliton. 

It is essential to mention some important character- 
istic features of the result shown in Figl?! Considering 
only the defect modes whose power is less than that of 
the fundamental mode (-Pa), one can find two regions 
along the axis Vq: regions of existence of the defect modes 
(white boxes) and regions where there are no solutions 
for P < Pa (gray boxes). The width of these regions 
can be controlled by changing the parameters of the fun- 
damental soliton as well as the width of the defect d. 
As discussed earlier, interaction of the high-amplitude 
soliton with the defect results in the reduction in width 
of the region of total reflection, as can be confirmed in 
FigEfb) where increasing the soliton amplitude results 
in the decrease in width of the corresponding regions of 
non-existence of solutions. 

Further analysis of the small-amplitude soliton re- 
vealed that the gray boxes exactly correspond to the in- 
tervals from the RCT diagram in Fig (31^ a) where the dis- 
crete soliton is completely reflected and the white boxes 
are in excellent agreement with regions where trapping 
of the discrete soliton by the defect occurs. In order 
to confirm this, we compare in FiglH] the capture coeffi- 
cient calculated numerically (blue line with dots) to the 
branches of existence of the OS and OA defect modes 
(lines) calculated as steady-state solutions of Eq. ([T5|). 
In order to prove that total transparency occurs at the 
matching points, in Fig llOl we compare the profiles of 
the defect modes, \u„\'^, taken at the points B, C and 
D (FiglS]) for ujs ~ —0.05 with snapshots of the soliton 
from direct numerical integration of the dynamical equa- 
tion ([Ij. The strength of the defect in the dynamical 



problem corresponds to that of the respective matching 
point and the time of the snapshot has been chosen to 
show the best symmetrical situation. In all these cases 
complete transmission of the soliton through the defect 
[similar to FiglSlJb)] has been observed. 
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FIG. 10. (Color online) Comparison of the profiles of the de- 
fect modes calculated from the stationary defect mode anal- 
ysis (red circles) and from direct numerical simulation (black 
squares) at the transmission points coinciding with the points 
of the maximum transmission (B, C, D). 



V. SPLITTING OF THE SOLITON BY THE 
POTENTIAL BARRIER: Vb < 

Let us now concentrate on the interaction of the dis- 
crete soliton with a potential barrier. Contrary to the 
previous case, now the soliton can only be either reflected 
or transmitted, depending on the relation between its ini- 
tial kinetic energy and the peak energy of the barrier. 
We determine the condition for the splitting of the ini- 
tial soliton into two identical parts after interaction with 
the impurity. To describe the dynamics of the discrete 
soliton and its interaction with the linear defect includ- 
ing the critical conditions for the splitting, we will use 
two conserved quantities of the DNLS equation which 
are the total power (number of particles) P determined 
in ([5]) and the Hamiltonian 



n I 7 , 



2 

n I J 



(16) 
(17) 

(18) 



which can be rewritten in terms of the frequency in the 
following form: 



W = — {Hq + Hd) = Wo + LUd. 



(19) 



Similar analysis was developed to study the scattering of 
bright solitons by a box-like barrier for the continuous 
NLS equation [ia |. 

In order to estimate the Hamiltonian (fT6|) for the single 
stationary soliton placed on top of the potential barrier 
taken in the form (fTTj) . we used the sech-type approx- 
imation "il^nit) = A sech(An)e*'^='* for the profile of the 
soliton. It should be stressed here that this limit is valid 
only if j4 1. In this limit the total power is simply 
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P = 2A. By substituting the sech-ansatz in 
resulting Hamiltonian (frequency) becomes: 



W the 



H = Ho{A) 



2AVof{A,d)i20) 



where the form of the function f{A,d) is determined by 
the width of the potential and the soliton amplitude and 
can be fitted with good accuracy using the following ex- 
pression 



f{A, d) = tanh{Ad + A/2 - O.i 



(21) 



having the two limits: as d — > the function / tends to 
A/ 2 and in the limit Ad > 1 the function / — > 1. In 
Fig llir a) we compare the function / calculated numeri- 
cally using the sech-profile ansatz with the fitting func- 
tion from (j2ip which shows good agreement for a large 
range of d. 

By fixing the width of the potential barrier and in- 
creasing its amplitude one can find the threshold Vcr 
when the discrete soliton breaks into two equal indepen- 
dent solutions moving from the center in opposite di- 
rections. The power of each part is P' = P/2 and the 
amplitude of the separated solitons is A' = A/2. Very 
close to the critical value the velocity of the split solitons 
is very small and its contribution to the total Hamilto- 
nian is negligible. By neglecting this contribution we can 
find the critical amplitude of the potential by compar- 
ing the Hamiltonian with the Hamiltonian calcu- 
lated for the pair of solitons far away from the defect: 
H' = 2Ha{A') = 4/3(A')3 = A^/&. Balance between 
these two Hamiltonians, namely, H = H' gives the criti- 
cal amplitude of the defect for which the splitting occurs: 



Vcr{A,d) 



A^ 



4/(A,d) 



(22) 



In FigfTTTb) we compare the critical value Vcr calculated 
analytically from with the corresponding value ob- 
tained from the dynamical problem. 

Below the critical amplitude of the potential, Vq < Vcr, 
one can observe breathing dynamics of the soliton placed 
on top of the defect. By exceeding the critical amplitude 
Vq > Vcr the soliton splits into two solitons and extra 
energy coming from the difference H — H' transforms into 
the kinetic energy (velocity) of the separated solitons. In 
Fig 112] these two cases are presented. 

Substituting back the critical amplitude ([22]) into 
and (fT9)) gives us the critical energy (frequency) at which 
splitting occurs: 



H. 



(cr) 



2 ' 



(cr) _ 



1 



= -A 



(23) 



Returning to the problem of the scattering of the soli- 
ton by the potential barrier, and considering now a mov- 
ing soliton placed far from the defect, we can again write 
its Hamiltonian which now has the form 
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FIG. 11. (Color online) In (a) comparison of the function / 
calculated analytically from (|2ip (red solid line) and numeri- 
cally as Hd/ {2AVo) (black dashed line with squares). In (b) 
comparison of the critical amplitude of the potential at split- 
ting as calculated analytically from (|22|l (red solid line) and 
from direct numerical integration of the dynamical equation 
(black dashed line with squares). The parameters used are 
A = 0.224, cj, = -0.05. 




FIG. 12. Density plot of the evolution of the soliton with 
ujs ~ —0.05 placed on top of the defect with width d = 2 
and amplitude (a) Vo = 0.027 < Vcr and (b) Vo = 0.028 > 
Vcr- The critical value of the potential amplitude for these 
parameters is Vcr ~ 0.028 from the analytical calculation and 
Vcr = 0.0273 from the direct numerical calculation. 



H — Hq + 



(24) 



where according to (HU]) ~ Puj^, = Av^ /2. Consider- 
ing the balance equation between the Hamiltonian of the 
moving soliton far from the defect ([2^ and the Hamilto- 
nian of the soliton placed on top of the defect at critical 
amplitude Her = H{Vo = Vcr), from ([TO]) we get the crit- 
ical velocity condition, Vcr = A, at which the separation 
between the reflection and transmission regimes occurs 
(due to the very large sensitivity to the initial condition 
it is very difficult to find numerically the point where the 
soliton splits into two identical parts). In Fig [13] the re- 
sult for the reflection (transmission) coefficient as a func- 
tion of (Vo,d) is shown. The initial soliton velocity is 
taken as the critical value v = Vcr and the separation 
between the reflection and transmission regions (bound- 
ary between the black and white regions) is in excellent 
agreement with the analytic prediction (dashed line) cal- 
culated from 
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R 



T 

■ -o;o4 ■ -o;o3 y -0:02 -O.'Ol 

FIG. 13. Density plot of the reflection coefflcient R and trans- 
mission coefficient T as a function of (Vb,d) with the initial 
soliton, taken at tOs — —0.05, placed at ns{t = 0) = —100 
and with initial velocity v = Vcr ~ A. The black dashed line 
corresponds to the region of splitting calculated from (P^. 



VI. CONCLUSION 

We have investigated the interaction of discrete soh- 
ton with a hnear potential taken in the form of a Gaus- 
sian function for both the smah-amphtude and high- 



amphtude regime where the sohtons are still mobile. In 
the case of an attractive defect, we observe the existence 
of windows of complete transmission or reflection as well 
as regions of partial capture of the soliton by the im- 
purity. The numerical findings are explained considering 
the stationary problem and the defect modes of the linear 
defect. Each peak of the total transmission corresponds 
to the matching of the power of the incoming discrete soli- 
ton with the power of the corresponding defect mode cal- 
culated from the underlying nonlinear eigenvalue prob- 
lem. Considering the interaction of the discrete soliton 
with a potential barrier we found analytically the criti- 
cal amplitude for the strength of the potential at which 
the splitting of the soliton into two identical solitons oc- 
cur. The numerical calculations show excellent agree- 
ment with the analytic prediction. We expect that the 
obtained results can be interesting for the construction 
of precise filters for solitons as well as for controllable 
switching devices in information processing in all-optical 
circuits. 
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